Enhanced diffusion due to active swimmers at a solid surface 
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We consider two systems of active swimmers moving close to a solid surface, one being a living 
population of wild-type E. coli and the other being an assembly of self-propelled Au-Pt rods. In 
both situations, we have identified two different types of motion at the surface and evaluated the 
fraction of the population that displayed ballistic trajectories (active swimmers) with respect to 
those showing random-like behavior. We studied the effect of this complex swimming activity on 
the diffusivity of passive tracers also present at the surface. We found that the tracer diffusivity is 
enhanced with respect to standard Brownian motion and increases linearly with the activity of the 
fluid, defined as the product of the fraction of active swimmers and their mean velocity. This result 
can be understood in terms of series of elementary encounters between the active swimmers and the 
tracers. 

PACS numbers: 87.16.-b, 05.65. -|-b 



Since the pioneering work of Wu and Libchaber [I] 
considerable efforts have been made to understand hy- 
drodynamic properties of active suspensions. Generally 
speaking, this is the name borne by fluids laden with self- 
swimming entities such as bacteria [2HS] , algae [7] or 
collections of active artificial swimmers [8] . Assemblies of 
microscopic motors dispersed in a fluid display emergent 
properties that differ strongly from passive suspensions. 
The momentum and energy transfer balances as well the 
constitutive transport properties are deeply modified by 
the momentum sources distributed in the bulk [21 [9]. 
Some of these anomalous properties have already been 
identified such as active diffusivity pQ [5] , anomalous vis- 
cous response [71 [TO], active transport and mixing [TT] 
as well as the possibility to use fluctuations to extract 
work |12j . The presence of living and apparently gregari- 
ous entities also offers the possibility to move collectively 
and organize at the mesoscopic or macroscopic level in 
the form of flocks and herds [HI [Il|. Similar collective 
effects were also identifled in suspensions of self-propeled 
inorganic particles |15| . In the bulk, swimming bacteria 
with flagella such as E. coli create in the far fleld-limit, 
a force-dipole velocity field and consequently, experience 
a hydrodynamic attraction toward surfaces !l6j. Then, 
it has been observed that E. coli smooth out their run- 
and-tumble movement and spend long times parallel to 
the surface undergoing circular motion as a consequence 
of the torque-free condition [171 [Hj ■ When the concen- 
tration becomes large, the E. coli population eventually 
associates collectively to form a bio-film. Even in the 
low concentration limit, the quantitative analysis of the 
near surface motion increase tremendously in complex- 
ity, a reason being the close field hydrodynamic forces 
that become prevalent and require a complex treatment 
of the lubrication hydrodynamic fields. However, even 
in this frame of description, it remains unclear whether 



the motion close to the surface is hydrodynamically sta- 
ble and if the presence of thermal noise is essential to 
account for the bacterium dwelling time at a surface (1 5]. 
Beyond the hydrodynamic interactions, more complex in- 
gredients may come into play such has the surface inter- 
action potentials (electrotastic or van der Waals) [2^ or 
more refined details of the bacterium physiology such as 
swimming speed variations and desynchronization dur- 
ing bacteria cell cycle [HI [22] . From the perspective of 
providing a fully consistent treatment of active hydrody- 
namics, with important applications for understanding 
bacterial transfer in biological micro- vessels, microfiuidic 
devices, or the formation of bio-films, a reliable descrip- 
tion of fiuid activity in the vicinity of a solid surface is 
strongly needed. In this letter, to tackle this open and 
timely question, we compared the behavior of two kinds 
of active micrometric swimmers: wild type E. coli K12 
and artificial self-propelling rods [8] , with completely dif- 
ferent propulsion mechanisms. In both cases, we moni- 
tor the swimmers' motions and their ability to activate, 
beyond Brownian motion, passive tracers, hence charac- 
terizing the active momentum transfer to the fiuid. 

Following the experimental procedure described in Ref. 
[16|, wild type E. coli K12 were grown overnight in rich 
medium (LB). After washing, they were transferred into 
MMAP, a motility medium supplemented with K-acetate 
(0.34 mM) and polyvinyl pyrohdone (PVP: 0.005%). 
They were incubated for at least an hour in that medium 
and, in some cases, so-called "baby cells" were selected by 
centrifugation and resuspended in MMAP. To avoid bac- 
terial sedimentation (isodense conditions), PercoU was 
mixed with MMAP (Ivol/lvol). We checked that under 
these conditions, the suspending fiuid was still Newtonian 
(viscosity 77=1.28x lO^^Fa s at 22°C). The overah con- 
centration of bacteria was controlled such as to prepare 
suspensions between 10^ and 10^" bact/ml. To study the 
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FIG. 1. Identification of the swimmer populations by tracking "active" swimmers (red tracks) and "random" swimmers (blue 
tracks), 0a is the corresponding fraction of active swimmers. Figs. l(a,b,g) correspond to E. coli (see inset in Ic) and Figs. 
l(d,e) to Au-Pt rods (see inset in If). The round black circles in (a,b) are 2/im latex beads; the white small circles in (d,e) are 
1/im Dynal beads. Figs, (c) and (f) display the relation between the active swimmers mean velocity Va and 0a. On Fig. 1(c), 
6 independent experiments with E. coli: IN cells (brown, red and green), mixture of IN and 2N cells (black) and 2N cells (pink, 
blue). Labels (1) and (2) are for A'^ = 100 and N = 200 cells in the observation field respectively. Fig. 1(f), for Au-Pt rods with 
varying proportions of inactive Au rods (2 independent experiments). Fig. 1(g), density probability of the observed bacterial 
tracks in the (TVc, {\0\)) space. The colormap goes from blue for vanishing probability to red for high probability. Two clusters 
are identified, centered at (0.9,0.3) and (0.17, 7r/2), corresponding respectively to the "active" and "random" swimmers. 



effect of bacterial activity on the diffusivity of passive 
tracers, latex beads of 1 or 2 diameter (Beckman- 
Coulter, density p=1.027 g/ml) were added to the sus- 
pensions. Experiments were performed in 110 fim thick 
chambers, built with two horizontal microscope cover 
slips separated by a glass spacer. To avoid sticking, cover 
slips were coated with PVP. The biological sample con- 
sisted in a drop of liquid (20 /il) placed between the two 
slides. The suspension was observed under an inverted 
microscope (Zeiss-Obzerver, Zl-magnification 40X) con- 
nected to a digital camera. The observation field AV^ was 
96x 128/im^ and 5 /ini in depth. In a first series of experi- 
ments, we measured the bacterial density profile through 
the entire height of the chamber. We obtained profiles 
similar to the ones published by Berke and coworkers [TB] , 
namely a flat density in the bulk and a strong density in- 
crease near the surfaces. However, the wild-type E. coli 
we used was significantly less attracted by the surfaces 
(2.5 times increase in density within 10 fim of the sur- 
face) than a mutant E. coli strain that does not display 
tumbling motion [T71 [TB] . Another series of experiments 
was performed with bimetallic Au-Pt self-propelled rods 
(length 1.2 /im and diameter 0.4 /im) that are very similar 
in size to the E. coli cell body (1 to 2 /im long, 0.8 /im di- 
ameter) but have no flagella (15 /im long for E. coli) (see 
inserts in Figs, [ij: and[lj'). They also have a much higher 
density {p—17 g/ml). In the presence of 1 to 10% hydro- 
gen peroxide, these particles are propelled in the axial 
direction towards the platinum end by the catalytic de- 



composition of the peroxide fuel [S] . Recent experiments 
and simulations are consistent with self-electrophoresis as 
the dominant propulsion mechanism |23l 124) . Here, the 
mode of propulsion is intrinsically different from the flag- 
ellar one. The experiments were conducted in a similar 
fashion to those involving E. coli, but in an open cham- 
ber (without the upper wall) , in order to allow the oxygen 
bubbles produced in the reaction to escape from the cell. 
The concentration of H2O2 was varied, as well as the 
concentration of active rods (n=3-20xl0^ rods/ml). We 
also used passive tracers (1/im diameter beads, Dynal- 
My-one, density p — 1.8 g/ml) or 2 /im diameter latex 
beads (Beckman-Coultcr, density p=1.02 7g/ml)) to fol- 
low the activation of the fluid by the Au-Pt rods. In all 
cases, all the particles in the suspension were localized at 
the bottom of the chamber due to sedimentation. Short 
videos (20s duration at 20/ps) were used to track both 
bacterial and self-propelled rods motion. 

In the following, we only focus on bacteria and rods 
moving close to the surface (less than 5pm). In both 
cases, we observed that not all swimmers display similar 
trajectories. This was expected as for wild type bacteria, 
the run or tumble dynamics may depend strongly on the 
microenvironment or on the position in the cell cycle. For 
Au-Pt rods, this is also consistent with previous observa- 
tions that even within a single batch, electrochemically 
grown rods have a range of catalytic activitylH]. We de- 
veloped a tracking program to analyze the short videos 
and obtained tracks for each swimmer present in the field. 
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We identified two major types of motion: a ballistic and a 
random one (see Fig.l(a,b,d,e)), and the swimmers that 
follow these motions are called "active" swimmers and 
"random" swimmers, respectively. To discriminate sys- 
tematically all the tracks, two parameters were defined. 
The first parameter {\9\) is the mean angle between two 
successive steps. For example, {\6\) — for straight tra- 
jectories and {\9\) = tt/2 for a purely random walk. The 
second discriminating parameter is based on the minimal 
circle diameter L that encompasses a given trajectory of 
duration T. For an acquisition time St (l/20s) and a 
mean step size Sr, the number Nc = is computed. 
When Nc is close to 1, the trajectory is associated with a 
straight line, whereas when Nc is small, its value points 
to diffusive motion. Therefore, each track is associated 
with these two numbers and in the {Nc, {\d\)) parameter- 
space, we could identify two clusters that clearly differ- 
entiated the active and random swimmers (see Figjljg)). 
Nevertheless, for very small or interrupted trajectories, 
the separation procedure remained ambiguous, so we sys- 
tematically discarded tracks shorter than 10 steps. In 
the case of artificial swimmers, we also managed to con- 
trol a priori the fraction of active swimmers by adding 
inactive rods (made only of gold) and keeping the to- 
tal number of rods at the surface constant. According to 
the trajectory classification, a fraction of active swim- 
mers was determined. Thus, for a mean number (N) of 
swimmers, identified in the field of vision, we define a 
density of "active swimmers" as ua = (f) a{N) / ISV . In 
Fig. [l] we display tracks during a time lag t — 1.5 s, for 
two populations of swimmers (a-b, bacteria and d-e, Au- 
Pt rods), having different <j)A (a,d small active- fraction 
and b,e high active fraction). In Figsjljc) and[T|f), we 
present the mean velocities of active swimmers Va as 
a function of their fraction (pA, for different experiments 
with bacteria and active rods. In the case of the bacterial 
suspension, we also tried several synchronization proto- 
cols to select bacteria, at different position in the cell 
cycle, showing different swimming characteristics. We 
were able to produce "baby-bacteria" populations (IN 
short cells: 1.12 /xm long) which were found to have a 4) a 
larger than the more mature bacteria populations (2N 
long cells: 2.5 /zm long). We took advantage of this dif- 
ference to look at the influence of (t)A on Va in bacterial 
suspensions. On Fig. 1(c) we display 6 independent ex- 
periments performed with populations having a majority 
of IN cells (brown, red and green), a mixture of IN and 
2N cells (black), a majority of 2N cells (pink, blue). For 
each sample, (pA was taken from suspensions showing an 
average of 100 or 200 bacteria in the observation field 
(represented by (1) and (2) on Fig. [Ij:, respectively). 0^ 
shows little influence on Va but Va could be different 
according to bacteria position within the cell cycle. In 
the case of the Au-Pt rods, Figjljf) shows a stronger in- 
dependence of Va on (pA, which is changed by varying 
the proportion of inactive Au-rods. This low dependence 
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FIG. 2. Enhanced diffusivity Dp of passive tracers measured 
as a function of J a ("activity flux"). Squares and circles sym- 
bols represent tracer of 2/im and l^m diameters, repsectively. 
Fig. 2(a) corresponds to the bacterial suspensions; IN cells 
(red and green) and unsynchronized (black, blue and pink). 
Each color defines an experiment performed over a range of 
bacterial dilution. Fig. 2(b) corresponds to suspensions of Au- 
Pt-rods: Mixture of active and inactive Au-rods (red circles, 
brown squares), and various H2O2 concentration from 2,5% to 
20%, (green, blue, black circles). The dashed lines are linear 
fits and the error bars, standard deviations. 



on (f)A is due to the low swimmer concentration of the 
suspension, where no collective behavior is observed. 

In the following, we will relate the motion of the pas- 
sive tracers to the number of active swimmers and their 
mean velocity. Passive tracer trajectories were analyzed 
(about 10 tracks in 300-image video at Ifps for bacteria 
and 40s sequences at 8fps for Au-Pt rods). No signif- 
icant stickiness between the spheres and the swimmers 
was observed and those rare cases were eliminated from 
the analysis. From these tracks, the mean passive tracer 
diffusion coefficient Dp was extracted consistently using 
two independent methods. The first one applied mean 
square displacements at long times (diffusive regime !.!]) 
of individual particles; the second one used particles as 
pairs in order to eliminate residual drift. In Fig. 2 the 
passive tracer diffusivities are displayed for all the exper- 
iments presented in Fig.l(c,f). Dp values are displayed 
as a function of J a — tiaVa, that we call the "activity 
ffux" . For experiments performed with the same tracer 
size, we observe a collapse of all data onto a linear curve: 



Dp = D\ 



13 Ja 



(1) 



where Dp is the Brownian diffusivity of the latex par- 
ticles in the vicinity of the surface, in the absence of 
swimmers. Note that due to lubrication forces, this value 
is smaller than the Brownian diffusivity expected in the 
bulk {Db = kBT/^TTrjd) and, for the parallel motion, 
they are related by Dp = aDg, where a < 1 is the par- 



allel drag correction factor p5H?7] . The a factor depends 
on the bead distance to the surface, vanishing at con- 
tact and going asymptotically to one at large distances. 
The beads are not at a fixed distance to the surface, but 
they are distributed according to the Boltzmann's fac- 
tor exp{—m*gz/kBT), where m* is the buoyant mass. 
Therefore a must be averaged with this factor. In the 



4 



active rods experiments we obtained a=0.7 both for the 
the buoyant d=l /xm and d=2 /Ltm latex spheres, value 
that agrees with the theoretical prediction given above 
(a=0.64 for (i=l/im and a=0.74 d=2 fim). In the ex- 
periments with bacteria, the suspended beads are almost 
isodense but they sediment anyway. The experimental 
fit gives a — 0.85. This value allows us to infer the 
density mismatch Ap — 0.008 g/ml, which is consis- 
tent with the experimental preparation. The collapse 
holds also for bacterial populations at different matura- 
tion stages (IN, 2N or unsynchronized mixtures). From 
dimensional analysis of expression (jlj , it can be seen that 
the prefactor /3 is a length to the fourth power. It varies 
from 5/im'''=(1.5/xm)* for bacteria, to 13/xm^=(1.9/Ltm)'* 
for active rods, but seems to be almost independent of 
the passive tracer size. Close to the plates the hydrody- 
namic perturbations created by the swimmers decay as 
the inverse cube of the distance, faster than in the bulk 
[28]. Therefore, at low concentrations, the enhanced dif- 
fusivity in Q, proportional to ua and Va, can be un- 
derstood as a result of a series of elementary encounters 
between active swimmers and the tracers: the number of 
encounters per unit time is proportional to uaVa- On 
the other hand, low Reynolds dynamics points out that 
the tracer displacement at each encounter is independent 
of the swimmer velocity and depends only on geometrical 
factors: the impact parameter, the swimmer dimensions 
and weakly on the tracer size through the Faxen correc- 
tion of passive transport ^29j, lSOj . The f3 factor comes out 
from averaging the tracer's displacements but its compu- 
tation is difficult because it requires a correct modeling 
of the near field interactions between the swimmer and 
the tracer, taking into account the detailled swimmer ge- 
ometry and the effect of the surface. 

In conclusion, we have characterized active momen- 
tum transfer close to solid surface for two active sus- 
pensions (wild-type bacteria and artificial self-propelled 
swimmers). The effect was measured using the diffusion 
enhancement of a passive tracer. In spite of the a pri- 
ori complexity of the hydrodynamics and essential dif- 
ferences in the propulsion modes, we demonstrated that 
the effect emerges quantitatively in a similar way. The 
resulting diffusion coefficient is the sum of the Brownian 
contribution near the wall and an active part, propor- 
tional to the product of the density of active swimmers 
and their mean velocity at the surface. The proportional- 
ity factor, scaling as the 4*'* power of a micron size length, 
encompasses the details of momentum transfer for each 
swimmer and is found to be weakly (if at all) sensitive 
to the probe diameter. Importantly, discriminating be- 
tween so-called "active" and "random" swimmer trajec- 
tories was crucial for predicting the induced transport 
phenomenon and we have developed a protocol to make 
such a distinction. The functional dependence of the en- 
hanced diffusivity is explained in terms of successive in- 
teractions between a single swimmer and the tracer, each 



one producing a net displacement. Our results justify to 
pursue a quantitative determination of such encounters 
based on simple hydrodynamic models [29] . 
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